Quantum spin nematics, dimerization, and deconfined criticality in quasi-one 

dimensional spin-1 magnets 
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We study theoretically the destruction of spin nematic order due to quantum fluctuations in quasi- 
one dimensional spin-1 magnets. If the nematic ordering is disordered by condensing disclinations 
then quantum Berry phase effects induce dimerization in the resulting paramagnet. We develop 
a theory for a Landau-forbidden second order transition between the spin nematic and dimerized 
states found in recent numerical calculations. Numerical tests of the theory are suggested. 

PACS numbers: 



Studies of low dimensional quantum magnetism pro- 
vide a good theoretical platform to develop intuition 
about the physics of strongly correlated many particle 
systems. In this paper we study various theoretical phe- 
nomena in spin S = 1 quantum magnets with SU (2) 
invariant nearest neighbor interactions. Specifically we 
focus on spin nematic order in such quantum magnets 
and its destruction by quantum fluctuations. 

A general Hamiltonian describing such spin S = 1 
quantum magnets takes the form 



(1) 



<tj> 



In real materials the ratio K/ J is probably small; how- 
ever it has been proposed that arbitrary values of Kj J 
can be engineered in ultra-cold atomic Bose gases with 
spin in optical lattices 1 . We will focus exclusively on a 
rectangular lattice where the couplings J, K on vertical 
bonds are a factor of A smaller than those on the horizon- 
tal bonds. This model was studied numerically recently 
(for K > 0) in an interesting paper by Harada et s&. 
In the isotropic limit A = 1, they found that there is a 
first order phase transition from a collinear Ncel state to 
a spin nematic state (along the line J = 0) with order 
parameter 
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even though there is no ordered moment < S >— 0. This 
spin nematic state corresponds to the development of a 
spontaneous hard axis anisotropy in the ground state. 
When A is decreased from 1 to make the lattice rect- 
angular quantum fluctuations are enhanced. The Neel 
and spin nematic phases then undergo quantum phase 
transitions to quantum paramagnets. Interestingly it is 
found that the spin nematic phase gives way to a dimer- 
ized quantum paramagnet where neighboring spin-1 mo- 
ments form strong singlets along every other bond in the 
horizontal direction. Further the quantum phase transi- 
tion itself appears to be second order in violation of naive 
expectations based on Landau theory but similar to the 



situations studied in Ref. (HQ} for other phase transitions 
in quantum magnets. 

In this paper we provide an understanding of these phe- 
nomena. First we provide general arguments relating the 
spontaneous dimerization with one route to killing spin 
nematic order by quantum fluctuations. When applied 
to one dimension our arguments explain the absence in 
numerical calculations^ of the featureless quantum dis- 
ordered spin nematic proposed by Chubukov 6 for spin-1 
chains. Further we show that a putative direct second or- 
der quantum phase transition between the spin nematic 
and dimerized phases is described by a continuum field 
theory with the action 



S = I d 3 x\ (<9 M - iAfj) D\ 2 + r\D\' 2 + (3) 
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Here D is a complex three component vector and is a 
non-compact U(l) gauge field. The nematic phase occurs 
when D condenses while D is gapped in the dimerized 
paramagnet. This theory is an anisotropic version of the 
non-compact CP 2 model (NCCP 2 ). The two compo- 
nent version - the anisotropic NCCP 1 model - describes 
Neel-VBS transitions of easy plane spin- 1/2 magnets on 
the square lattice^. A second order nematic-dimer tran- 
sition on the rectangular lattice is possible if doubled in- 
stantons in the gauge field are irrelevant at the critical 
fixed point of this theory. These instantons are relevant 
at the paramagnetic fixed point of Eqn. [3l This leads to 
confinement of the D fields and to dimer order. Indeed 
the dimer order parameter is simply the single instanton 
operator. A direct second order nematic-dimer transi- 
tion is thus accompanied by the dangerous irrelevance 
of doubled instantons and the associated two diverging 
length/ time scales. 

We will provide two different arguments to justify our 
results. First we address the quantum disordering of the 
nematic based on general effective field theory considera- 
tions that focus on the properties of topological defects of 
the nematic order parameter. Second we provide a more 
microscopic argument based on an exact slave 'triplon' 
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representation 7 of the S = 1 operators. 

The order parameter manifold for the spin nematic 
state may be taken to be the possible orientations of 
the spontaneous hard axis d (the "director") and thus 
is S 2 /Z2- The Z2 simply reflects the fact that d and — d 
are the same state. The spin nematic state allows for 
Z 2 point vortex defects ( "disclinations" ) in two space di- 
mensions. The director d winds by tt on encircling such a 
disclination. As discussed by Lammert et al£ for classical 
nematics the fate of the disclinations crucially determines 
the nature of the "isotropic" phase obtained when the ne- 
matic is disordered by fluctuations. If the transition out 
of the nematic occurs without condensing the disclina- 
tions then a novel topologically ordered phase - inter- 
preted in the present context as a quantum spin liquid 
- obtains. However the nematic may also be disordered 
in a more conventional way by condensing the disclina- 
tions. In the present context we argue that non-trivial 
quantum Berry phases associated with the disclinations 
lead to broken translational symmetry in this quantum 
paramagnet. Furthermore this transition may be second 
order as described below (unlike the classical nematic- 
isotropic transition). 

Following Lammert et al^ we consider the quantum 
phase transition out of the nematic using an effective 
model in terms of the director d. The d, — d identifica- 
tion requires that the d vector is coupled to a Z2 gauge 
field. Thus we consider the following action on a three 
dimensional spacetime lattice: 

S = S d + S B (4) 
S d = - ^2 t^a^dr.dr+f, (5) 

<r,r+fj,> 

Here r represent the sites of a cubic spacetime lattice, 
/i = (x,y,T), cr M (r) = ±1 is a Z2 gauge field on the link 
between r and r + /i. The term Sb is the Berry phase to 
be elaborated below. 

The Berry phases arise from the nontrivial quantum 
dynamics of the d vector and can be understood very 
simply by considering a single quantum spin S = 1 with 
a time varying hard axis d(r) that represents the fluctu- 
ating local director field: 

H=(d{r).s) 2 (6) 

For a time independent d, the ground state is simply the 
state where the projection S.d = 0. The Berry phase is 
obtained by considering a slow time varying closed path 
of d in the adiabatic approximation. There are two kinds 
of such closed paths that are topologically distinct. First 
there are paths for which d returns to itself. For such 
paths it is easy to see that the Berry phase factor is 1. 
Then there are closed paths where d returns to — d. In 
the adiabatic approximation with S = 1 it is easy to see 
that the wavefunction acquires a phase of tt for such a 



path. Thus there is a Berry phase of —1 for closed paths 
where d returns to —d. 

The phase factor of —1 for nontrivial closed time evolu- 
tions of d at a spatial site may be naturally incorporated 
into the effective lattice model of Eqn. [5] above. First we 
note that a closed loop in time where d winds by tt corre- 
sponds to a configuration with Z 2 gauge flux —1 through 
the loop. The Berry phase is thus simply 

e~ s * =]Ja T (r) (7) 

r 

At each space point the product over the time-like bonds 
measures the flux of the Z 2 gauge field through the closed 
time loop at that point. Precisely this Berry phase factor 
arises in Z 2 gauge theoretic formulations of a number of 
different strong correlation problems 9 , and the theory is 
known as the odd Z2 gauge theory. Thus an appropriate 
effective model for disordering the S — 1 spin nematic 
state is a theory of d coupled to an odd Z2 gauge theory. 

The spin nematic ordered phase corresponds to a con- 
densate of d. In this phase the Z2 disclinations are simply 
associated with Z2 flux configurations of the gauge field. 
Thus the Berry phase term associated with the gauge 
field directly affects the dynamics of the disclinations. 
Disordered phases where d has short ranged correlations 
may be discussed by integrating out the d field. The re- 
sult is pure odd Z2 gauge theory on a spatial lattice with 
rectangular symmetry. This theory is well understood. 
It is conveniently analysed by a duality transformation 
to a stacked fully frustrated Ising models followed by a 
soft spin Landau-Ginzburg analysis^. This leads to a 
mapping to an XY model with four-fold anisotropy: 

S v = -t v ^2 cos(<j> R — <j>Ri)— cos{A4r) (8) 

<RR'> R 

Here R, R' are sites of the dual cubic lattice. The 
real and imaginary parts of the field e lc ^ R correspond to 
Fourier components of the Z 2 vortex near two different 
wavevectors at which the quadratic part of the Landau- 
Ginzburg action has minima. The anisotropy is four- 
fold on the rectangular spatial lattice as opposed to the 
8-fold anisotropy that obtains with square symmetry^. 
There is a disordered phase where the Z 2 vortex has 
short-ranged correlations: this corresponds to the topo- 
logically ordered quantum spin liquid in the original spin 
model. In addition there are ordered phases associated 
with condensation of the e % ^ R . These phases break trans- 
lation symmetry. For the rectangular lattice of interest 
the natural symmetry breaking pattern is dimerization 
along the chain direction. 

We thus see that Berry phases associated with the 
quantum dynamics of the director d lead to dimeriza- 
tion when the nematic order is disordered by condensing 
the Z2 disclinations. This analysis can be easily repeated 
in one spatial dimension. Then the Z2 disclinations are 
point defects in spacetime. These are described by the 
odd Z2 gauge theory in 1 + 1 dimensions which is always 
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confined and which has a dimerized ground stated. In 
particular this argument shows that for S = 1 chains a 
featureless disordered spin nematic state will not exist. 

Returning to two dimensions we may now write down a 
field theory for the nematic-dimer transition. The Berry 
phases on the disclinations are encapsulated in Eqn. [8] 
We now need to couple these back to the d vector. The 
main interaction between d and e l ^ R is the long ranged 
statistical one: on going around a particle created by e 1 ^ 
the vector d acquires a minus sign. 

We proceed by first ignoring the k term in Eqn. [5] and 
using a Villain representation to let 

Sv= E u $uv (9) 

<RR'> 

The 2rr' are integer valued currents of the e 1 ^ that sat- 
isfy the conditions 

Vj = (10) 

(-iy = ]Ja (11) 
p 

The first equation expresses current conservation (which 
holds at k = 0). In the second the symbol FJ p refers to a 
product over the four bonds of the direct lattice pierced 
by < RB! >. This term ensures that an e 1 ^ particle acts 
as 7r flux for the d field. Solving Eqn. [10] by j = V X A 
(with integer) we get 

(-1)VXA = -Q a (12) 
p 

Now write A = 2a + s with a an integer and s = 0, 1 so 
that 

IB- 1 )" = 11* ( 13 ) 
p p 

which can be solved by choosing 

(-l) s = l-2s =a (14) 

The integer constraint on A may be implemented softly 
by including a term 

— tgcos(2ira) = —tga rr icos(irA rr i) (15) 

We now separate out the longitudinal part of A by letting 

A^A+-V6 (16) 

7T 



After a further rescaling A 

S = 
S e = 

S A = 



^ we finally get the action 

(17) 
(18) 

(19) 



Sd + Sg + Sa 

—to a rr 'Cos(6 r — 6 r > + A rr i) 



U 



with Sd given in eqn. [4] The sum over a can now be per- 
formed. The universal prperties are correctly captured in 
a model that keeps the lowest order cross term between 
t d and te. We therefore get 

S = - t^cos(9 r - r i + A rr ,)d r .d r i + U Y^(V x A) 2 

<rr'> P 

(20) 

with ~ fjjto. It is instructive to introduce the complex 
vector D r = e r d r that satisfies 

l^| 2 = 1 (21) 
D x D* = (22) 

The second condition may be imposed softly by including 
a term 

v\DxD*\ 2 = -v({Df(D*) 2 -l) (23) 
with v > 0. Thus we arrive at the model 

S = S d + Sa (24) 
S D = -t ^ e lA ^' 3;.D,, + c.c- v(D) 2 (D*) 2 (25) 

<rr'> 

with D satisfying |75| 2 = 1. Atw = this is the lattice 
CP 2 model with a global SU(3) symmetry associated 
with rotations of the D field. The v term breaks the 
SU(3) symmetry down to SO(3). Eqn. [3] is precisely a 
soft-spin continuum version of the lattice action above. 

We now consider the role of the four fold anisotropy on 
the disclination field e 1 ^ (the k term in Eqn [8]). Without 
this term the number conjugate to </> is conserved. In 
the dual description this translates into conservation of 
the magnetic flux of the U(l) gauge field A. Thus at 
k = the gauge field is noncompact. The n term however 
destroys this conservation law - indeed four disclinations 
can be created or destroyed together. In the effective 
model of Eqn. GH a disclination in D corresponds to a 
configuration where the gauge flux is equal to tt. Thus 
the k term may be interpreted as a doubled 'instanton' 
operator that changes the gauge flux by 47r. 

The nematic order parameter is simply related to the 
D fields: 

Q „ = (£M±±i±.^ (26) 

Thus when D condenses nematic order develops. The 
paramagnetic phase occurs when D is gapped. In the 
absence of instantons the low energy theory of this phase 
has a free propagating massless photon. Instantons how- 
ever gap out the photon and confine the D fields. The 
dimer order parameter e 1 ^ is the single instanton operator 
and gets pinned in this phase. A direct second order tran- 
sition between the nematic and dimerized states can thus 
occur if doubled instantons are irrelevant at the critical 
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fixed point of the anisotropic NCCP 2 action associated 
with the condensation of D. 

A different more microscopic argument can also be 
used to justify Eqn. [3] and provides further insight. Con- 
sider the following exact representation^ of a spin-1 op- 
erator at a site i in terms of a 'slave' triplon operator 

Wi. 



Si = 



iw\ x Wi 



(27) 



together with the constraint wj.Wi = 1. The Wi satisfy 
usual boson commutation relations. The nematic order 
parameter is readily seen to simply be 



ia/3 



S a /3 _ wjw/3 + C.C 

3 2 



(28) 



As with other slave particles, this representation leads to 
a U(l) gauge redundancy associated with letting 

Wi e la "Wi (29) 

at each lattice site. It is convenient to first consider the 
special point Jij = where the Hamiltonian in Eqn. [1] is 
known to have extra SU (3) symmetry. Then H may be 
rewritten (upto an overall additive constant) 



(30) 



This is invariant under a global multiplication of Wi by 
an SU (3) matrix U on one sublattice and by U* on the 
other. Such magnets were studied in detail in Ref. [HI and 
we can take over many of their results. A standard mean 
field approximation with < Wi.Wj >^ yields a param- 
agnetic phase with gapped w particles in the d = 1 limit 
while in two dimensions the w condense thereby breaking 
the St/ (3) symmetry. The theory of fluctuations beyond 
mean field includes a compact 1/(1) gauge field. In the 
paramagnetic phase instanton fluctuations of this gauge 
field confine the w particles and their Berry phases lead 
to dimerization on the rectangular lattice. The results 
of Ref. |12| now imply that the transition associated with 
w condensation is described by an NCCP 2 model with 
doubled instantons, i.e it is precisely of the form of Eqn. 
|3] but with v = 0. The triplon wJj on the A sublattice 
~ D while on the other sublattice Wi ~ D*. Thus we see 
that the Neel vector N is simply related to D through 



N 



-iD* x D 



(31) 



For the SU (3) symmetric Hamiltonian all eight compo- 
nents of the tensor D^Dp — \D\ 2 5 a p/3 have the same 



correlators. The symmetric part of this tensor is the ne- 
matic order parameter and the antisymmetric part is the 
Neel vector. 

If now a small J < is turned on the SU (3) symmetry 
is explicitly broken down to SO(3). This sign of J disfa- 
vors Neel ordering so that nematic ordering wins in the 
two dimensional limit. The NCCP 2 field theory of the 
transition to the dimer state must then be supplemented 
with an anisotropy term v\N\ 2 with v > which due to 
Eqn. [21] is precisely the anisotropy term of Eqn. [3] 

What may we say about the NCCP 2 field theory 
and the fate of doubled instantons? First as the criti- 
cal theory is relativistic the dynamical critical exponent 
2 = 1. Next we note that the instanton scaling dimen- 
sion is expected to be bigger for NCCP 2 as compared to 
NCCP 1 . In the isotropic case existing estimates^ give 
0.63 for the single instanton scaling dimension. In a 
naive RPA treatment of the gauge fluctuations the in- 
stanton scaling dimension scales like m 2 N where m is 
the instanton charge and N is the number of boson com- 
ponents. Thus within this approximation we estimate 
the doubled instanton scaling dimension in NCCP 2 as 
|(2) 2 (.63) « 3.78. This admittedly crude estimate nev- 
ertheless suggests that doubled instantons may be irrel- 
evant for NCCP 2 . 

The possible irrelevance of the doubled instantons 
has dramatic consequences for the phenomena at the 
nematic-dimer transition. It implies that the critical 
fixed point has enlarged t/(l) symmetry associated with 
conservation of the gauge flux exactly like in Ref. H. 
This enlarged symmetry implies that the (w, 0) columnar 
dimer order parameter may be rotated into the (0, 7f) 
columnar dimer or into plaquette order parameters at 
(0,7r), (tt, 0). Thus right at the critical point all these 
different VBS orders will have the same power law cor- 
relations. It will be an interesting check of the theory of 
this paper to look for this in future numerical calcula- 
tions. 

In summary we have studied the destruction of spin 
nematic order by quantum fluctuations in quasi-one di- 
mensional spin-1 magnets. We showed that Berry phases 
associated with disclinations lead to dimerization if the 
nematic is disordered by their condensation. We pre- 
sented a continuum field theory for a Landau-forbidden 
second order transition between nematic and dimerized 
phases that generalizes earlier work on other transitions. 
Future numerical work or cold atoms experiments may 
be able to explore the physics described in this paper. 
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